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ABSTRACT 

We present Hi line profiles for various models of circumstellar shells around red giants. 
In the calculations we take into account the effect of the background at 21 cm, and show 
that in some circumstances it may have an important effect on the shape and intensity 
of the observed line profiles. We show that self-absorption should also be considered 
depending on the mass loss rate and the temperature reached by circumstellar gas. 

Hi emission from circumstellar shells has been mostly reported from stars with 
mass loss rates around 10 -7 M 0 yr -1 . We discuss the possible reasons for the non 
detection of many sources with larger mass loss rates that are hallmarks of the end of 
the AGB phase. Although radiative transfer effects may weaken the line emission, they 
cannot alone account for this effect. Therefore, it seems likely that molecular hydrogen, 
rather than atomic hydrogen, dominates the composition of matter expelled by stars 
at the end of their evolution on the Asymptotic Giant Branch. However sensitive 
HI observations can still yield important information on the kinematics and physical 
properties of the circumstellar material at large distances from central stars with heavy 
mass loss, despite the low abundance of atomic hydrogen. 

Key words: stars: AGB and post-AGB - circumstellar matter - stars: individual: 
Y CVn - stars: mass loss - radio lines: stars. 


1 INTRODUCTION 

Asymptotic giant branch (AGB) stars are undergoing mass 
loss via slow stellar winds (Olofsson 1999). They get sur¬ 
rounded by expanding circumstellar shells that are well 
characterized by CO emission lines in the millimeter range 
(Knapp & Morris 1983, Loup et al. 1993). From the model¬ 
ing of the CO line profiles, expansion velocities are derived 
to be in the range from a few to 40 kms -1 , and mass loss 
rates, from ~ 10 -7 to a few 10 -4 M 0 yr -1 . 

Most of the matter in circumstellar shells should be 
in hydrogen (~70% in mass). Early searches for 21-cm 
line emission from atomic hydrogen were mostly unsuccess¬ 
ful (Zuckerman et al. 1980, Knapp & Bowers 1983) lead¬ 
ing to the conclusion that hydrogen in stellar envelopes is 
molecular rather than atomic. However, more recently, the 
Hi line at 21 cm has been detected from several AGB stars 
(Gerard & Le Bertre 2006). A number of the detected stars 
have also been imaged in HI with the Very Large Array 


(VLA, Matthews & Reid 2007, Matthews et al. 2008, 2011, 
2013). In general, the detected sources have mass loss rates 
~ lO -7 M 0 yr -1 , and they show a narrow Hi line profile 
(3-5 kms -1 ). On the other hand, many sources with large 
mass loss rate (Js lO -6 M 0 yr -1 ) remain undetected. 

Again, a possible explanation could be that hydrogen 
in sources with large mass loss rates is mostly in molec¬ 
ular form. However, this is not entirely satisfactory for at 
least two reasons. First, molecular hydrogen is expected to 
be photo-dissociated by UV photons from the interstellar 
radiation field (ISRF) in the external parts of circumstellar 
shells (Morris & Jura 1983). Second, the mass loss rate from 
AGB stars is supposed to be highly variable (e.g. Vassiliadis 
& Wood 1993), so that a given star should be surrounded by 
interacting shells, originating from different episodes of mass 
loss at different rates, from 10 -8 M 0 yr -1 upwards, and dif¬ 
ferent velocities (Villaver et al. 2002), with probably a mix 
of atomic and molecular hydrogen. 

Another possibility is that radiative transfer effects are 
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weakening the 21-cm line emission, thus rendering difficult 
the detection of H I. The narrowness of the line profiles could 
increase the optical depth at particular frequencies, thus ex¬ 
acerbating self-absorption. The presence of an intense back¬ 
ground at 21 cm may also play a role (Levinson & Brown 
1980). In this paper, we examine the formation of the HI line 
in circumstellar shells around AGB stars, taking into ac¬ 
count these effects. 


2 CONTEXT 

Hydrodynamical models of dust-forming long-period vari¬ 
able stars reproduce outflows with mass loss rate ranging 
from ~ 10 -8 to 10~ 4 M@ yr _1 , and expansion velocity rang¬ 
ing from a few to 40 kms' 1 (Winters et al. 2000). Close to 
the central star, the density is much larger than in the inter¬ 
stellar medium (ISM), and the wind expands freely (freely 
expanding wind). We hereafter refer to the case where the 
wind expands freely as scenario 1. However, at some point, 
the stellar wind is slowed down by the ambiant medium, and 
a structure develops in which both circumstellar and inter¬ 
stellar matter accumulates. Such structures were discovered 
with IRAS owing to the presence of spatially resolved 60/rm 
continuum emission (Young et al. 1993) and are referred to 
as ‘detached shellsQ- 

Subsequently, the Hi line at 21 cm was detected in 
the direction of several resolved IRAS sources (Le Bertre 
& Gerard 2004; Gerard & Le Bertre 2006). The Hi emission 
is more extended than the CO emission and generally corre¬ 
sponds with the far-infrared size. Also the main components 
of the Hi line are narrower than the CO line profiles, thus 
demonstrating that stellar outflows are slowed down by an 
external medium, and that the Hi line profiles can be used 
to probe the kinematics in the external regions of circum¬ 
stellar shells. Imaging obtained with the VLA has revealed 
detached Hi shells with narrow line profiles, sometimes as¬ 
sociated with extended tails (Matthews et al. 2013). 

Libert et al. (2007) developed a model of a detached 
shell (hereafter scenario 2) adapted to the HI emission from 
the prototypical case of Y CVn. In this model, a double¬ 
shock structure develops, with a termination shock, fac¬ 
ing the central star, where the supersonic stellar wind is 
abruptly slowed down, and a leading shock (i.e. outward 
facing shock) where interstellar matter is compressed and 
integrated in the detached shell. These two limits define the 
detached shell which is thus composed of circumstellar and 
interstellar matter separated by a contact discontinuity. The 
circumstellar matter is decelerated, by a factor ~ 4, when 
it crosses the termination shock, and heated. It then cools 
down, while the expansion velocity further decreases and 
the density increases until the contact discontinuity. The 
model assumes spherical symmetry and stationarity. Adopt¬ 
ing an arbitrary temperature profile, the equation of motion 
is solved numerically between the two limits of the detached 
shell. The Hi line profile is then calculated assuming that 
the emission is optically thin. 


1 This designation has also been used for shells resulting from an 
outburst of mass loss and observed in CO molecular lines (Olofs- 
son et al. 2000). In this paper, we adopt to the definition of Young 
et al. 


Libert et al. (2007) could reproduce the Hi line profile 
observed in Y CVn and generalize the model to other cases 
of sources with mass loss rates of about 10~' Mgyr -1 (e.g. 
Libert et al. 2010). However, as discussed in Matthews et 
al. (2013), although this model is effective in reproduc¬ 
ing spatially integrated line profiles, it faces some difficul¬ 
ties in reproducing the spatial distribution of the Hi emis¬ 
sion. Namely, at small radii, it seems to predict too much 
flux in both the shell component and the freely expanding 
wind component of the spectral profile. The mass loss vari¬ 
ations with time and the interaction between consecutive 
shells may need to be taken into account (see Villaver et al. 
2002). The more elaborate numerical modeling performed 
by Villaver et al. follows the mass loss and velocity varia¬ 
tions of the stellar winds of thermally pulsing AGB stars, 
and describes self-consistently the formation and develop¬ 
ment of circumstellar shells embedded in the ISM (hereafter 
scenario 3). They show that large regions (up to 2.5 pc) of 
neutral gas may be formed around AGB stars. 

In this work, we will simulate the Hi line profiles that 
would be expected for these three kinds of scenarios and ex¬ 
amine how they compare with observations presently avail¬ 
able. 


3 GENERAL COMMENTS ON THE 

FORMATION OF THE H i LINE PROFILE 


The Hi line at 21 cm is produced by a transition between 
two levels of the ground state of hydrogen. The rest fre¬ 
quency, vo, corresponds to 1420.4 MHz, and the transition 
probability, A 21 , to 2.87 10 -15 s^ 1 . As the transition proba¬ 
bility between the two levels is very low, the natural width of 
the line is extremely small, and the frequency dependence of 
the absorption and emission coefficients are entirely defined 
by the atomic hydrogen velocity. This is determined by the 
bulk motion of the gas (including, if present, turbulence), 
and the hydrogen kinetic temperature (T). We assume local 
thermal equilibrium (LTE). Also, in general, as the frequen¬ 
cies are small, one may replace exp(hv/kT) — 1, by hv/kT 
( hv/k ~ 0.07K). 

The intensity emitted, at a frequency v , along an ele¬ 
ment ds, dl v , is given by : 

dlv hv f \ a / -1 \ 

-dT = i ^ 2 ^ 21 (1) 

where n 2 (v) is the number density of hydrogen atoms in 
the upper level with a velocity along the line of sight corre¬ 
sponding to v. 

Similarly the absorption part is given by : 

t < 2 > 

with the absorption coefficient, k{v), given by : 


k(v) 


c?n\(v)g 2 . hv 
8nvo 2 gi 21 kT 


( 3 ) 


(cf. equation 3.14 in Lequeux 2005). (?2 and gi are the statis¬ 
tical weights of the upper and lower levels, respectively, and 
ni ( v ) is the number density of hydrogen atoms in the lower 
level with a velocity along the line of sight corresponding to 
v. It is interesting to note that, in general, the absorption 
coefficient varies as 1/T. 
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For the Hi line at 21 cm: g\ = 1, g 2 == 3, and for a 
population in statistical equilibrium, m = 1/4 uh, and n -2 = 
3/4 tih, where nn is the total number density of hydrogen 
atoms in the ground state ( uh = ni + 122 ), and assuming 
exp(hu/kT) = 1 for the Boltzmann factor. 

In practice, spectra are represented as a function of the 
velocity, V. We may thus replace ni(u) by l/4riH(V) and 
n 2 [y) by 3/4 uh{V). With this convention, the absorption 
coefficient can be rewritten as: 


n{V) 


3(?uh{V) . h 
32-rvq 21 kf' 


(4) 


Finally the radiative transfer equation is written as : 


dI(V) 

ds 


3hu 0 

167T 


n H (V)A 2 1 - 


3 c 2 n H (V) 
32m/o 


A ^nn 


(5) 


At radio frequencies, it is usual to express the intensity 
in terms of the equivalent temperature of a blackbody that 
would give the same intensity in the same spectral domain. 
With this convention, the boundary condition can be defined 
through a background brightness temperature, Tbg ■ 




Figure 1. Density and temperature profiles for an out¬ 
flow in uniform expansion (scenario 1, V eX p = 10 km s' 1 , 
M = 10“ 5 M 0 yr —1 ). 


I + (V) 



Tbg(V), 


( 6 ) 


4 SIMULATIONS 


I + refering to the incoming intensity on the rear side of the 
shell. The background is the sum of the 3-K cosmic emis¬ 
sion, the synchrotron emission from the Galaxy and the 
Hi emission from the ISM located beyond the circumstel- 
lar shell with respect to the observer. The first component 
is a continuum emission, which is smooth, angularly and 
spectrally. The second component is also smooth spectrally, 
but it presents a strong dependence with galactic latitude, 
and shows also some sub-structures. The sum of these two 
components has been mapped with a spatial resolution of 
0.6° by Reich (1982), Reich & Reich (1986) and Reich et al. 
(2001). The third component (Hi emission from the ISM) 
shows both strong spatial and spectral dependences, which 
make it a serious source of confusion. It has been mapped 
with a spatial resolution of 0.6° and a spectral resolution 
of 1.3 km s' 1 by Kalberla et al. (2005; Leiden-Argentina- 
Bonn, LAB, survey). Surveys of selected regions of the sky, 
in particular along the Galactic Plane, have been obtained 
at a resolution down to 1', and show spatial structures, like 
filaments or clouds, at all sizes (e.g. Stil et al. 2006). 

Away from the Galactic Plane, typical values range 
from Tbg ~3-5 K, outside the range of interstellar Hi emis¬ 
sion, to Tbg ~10-20 K inside an interstellar Hi emission. 
Close to the Galactic Plane, the continuum may reach 
Tbg ^10-20 K, and, including the interstellar Hi emission, 
the background may reach Tbg ~100 K. Note that the back¬ 
ground temperature, Tbg, is not directly related to the ki¬ 
netic temperature of the surrounding ISM. 

In addition, in some cases, a radio-source may be seen in 
the direction of a circumstellar shell. In such a case, we have 
an unresolved continuum emission (see e.g. Matthews et al. 
2008). Such a source may be useful to probe the physical 
conditions within the circumstellar shell in a pencil-beam 
mode. 


For the present work, we adapted the code developed by 
Hoai et al. (2014). It is a ray-tracing code that takes into 
account absorption and emission in the line profile. It can 
handle any kind of geometry, but for the purpose of this 
paper we restricted our simulations to circumstellar shells 
with a spherical geometry as described in Sect .[2] We assume 
that, in each cell, the gas is in equilibrium and that the 
distribution of the velocities is maxwellian. 

In this section, we explore the line profiles for a source 
that is not resolved spatially by the telescope, and assume a 
uniform response in the telescope beam (boxcar response, cf. 
Gardan et al. 2006). We also assume that the line profiles 
can be extracted from position-switched observations, i.e. 
that there is no spatial variation of the background. The 
flux densities are expressed in the units of Janskys (Jy), 
where 1 Jy=10' 26 Wm' 2 Hz' 1 . 

We performed various tests in order to evaluate the ac¬ 
curacy of the simulations. It depends mainly on the mass loss 
rate of the central source and on the size of the geometrical 
steps adopted in the calculations. For the results presented 
in this section, the relative error on the line profile ranges 
from ~ 10~ 6 , for mass loss rates of 10'' Mgyr -1 , to a few 
10“ 3 , for mass loss rates of 10' 4 Mgyr' 1 . 

4.1 Freely expanding wind (scenario 1) 

We consider a spherical wind in free expansion at 
Vexp=10kms“ 1 . The distance is set at 200 pc, and the mass 
loss rate is varied from 10' 7 to 10~ 4 Mgyr' 1 . We assume 
that the gas is composed, in number, of 90% atomic hydro¬ 
gen, and 10% 4 He. We assume a temperature dependence 
proportional to r' 0,7 , r being the distance to the central 
star, out to the external boundary (0.17 pc) where the tem¬ 
perature drops to 5 K. This temperature of 5 K is proba- 
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Figure 2. Hi line profiles of shells in free expansion for vari¬ 
ous mass loss rates with no background. The profiles for 10 7 , 
10 -6 , and 10~ 5 Moyr _1 are scaled by factors 1000., 100., and 
10., respectively. The distance is set at 200 pc. 


bly underestimated as the photoelectric heating by grains 
absorbing UV photons is expected to raise the tempera¬ 
ture of the gas in the cool external layers of shells around 
stars with high mass loss rate (Schoier & Olofsson 2001). 
On the other hand, temperatures as low as 2.8 K have been 
reported in some high mass loss rate sources (e.g. U Cam, 
Sahai 1990). Such low temperatures are only expected in 
the freely expanding regions of the circumstellar shells. The 
density and temperature profiles are illustrated in Fig. [l] for 
the lO - 5 M 0 yr -1 case. 

In a first set of simulations (Fig.0 we calculate the in¬ 
tegrated emission (within a diameter, <j> = 6') with no back¬ 
ground, in order to estimate the effect of self-absorption (in 
fact the background should have a minimum brightness tem¬ 
perature of 3K, cf. Sect. 0. Self-absorption starts to play 
a clear role for 10 ~ 6 Mgyr -1 , with an intensity that is re¬ 
duced, and a profile that changes its shape from almost rect¬ 
angular to parabolic. A slight asymmetry of the line profile 
is also present (although not discernible by eye in the figure), 
with more absorption on the blue side, due to the outwardly 
decreasing temperature, an effect which has already been 
described in the case of molecular emission from expanding 
circumstellar envelopes (Morris et al. 1985). 

In a second set of simulations, the mass loss rate is kept 
at 1CP 5 Mq yr -1 , and the background is varied from OK (as 
above), to 10K (Fig.0. The effects noted previously are am¬ 
plified by the background, in particular with an absorption 
developing on the blue side of the profile, and then extend¬ 
ing to the complete spectral domain when the background 
temperature reaches 10 K. 

We adopted a temperature dependence proportional to 
r ~ 0,7 which fits the results obtained with a radiative trans¬ 
fer model by Schoier & Olofsson (2001). A shallower depen¬ 
dence would increase the temperature in the outer layers of 
the circumstellar shell and thus reduce the effects of self¬ 


Figure 3. Hi line profiles of a shell in free expansion for 
M=10 -5 Mq yr —1 , and for various background levels ( Tbg= 0, 
3, 5, 7, 10K). 

absorption, as well as the absorption of the background ra¬ 
diation. An external source of heating (e.g. by photoelectric 
heating) would have the same influence. 

4.2 Single detached shell (scenario 2) 

We adopt the model developed by Libert et al. (2007). 
It has been shown to provide good spectral fits of the 
Hi observations obtained on sources with mass loss rates 
~ 10 -7 Mq yr -1 (cf. Sect.0. 

As in Sect. 14. ll we assume a spatially unresolved source 
at 200 pc with a mass loss rate of 10 -7 Mq yr -1 . The internal 
radius of the detached shell is set at 2.5' (or 0.15 pc). Simi¬ 
larly, we examine the dependence of the line profile for mod¬ 
els with various masses in the detached shell (M_dt,cs) and 
various background levels (Figs.0and0. The parameters of 
the four cases illustrated in Fig. [4] are given in Table 0 The 
free-wind expansion velocity is taken to be V exp = 8 kms~ 1 . 
At the termination shock the downstream temperature is 
given by T/ ~ (3 pm.sr)/(16 k) V exp 2 ~ 1800 I< (equ. 6.58 
in Dyson & Williams 1997) with rriH, the mass of the hy¬ 
drogen atom, and p, the mean molecular weight. For the 
temperature profile inside the detached shell we use the ex¬ 
pression (9) in Libert et al. (2007) with a temperature in¬ 
dex, a = —6.0. The temperature is thus decreasing from 
~ 1800K, to T/, at the interface with ISM, r/. The density, 
velocity and temperature profiles are illustrated in Fig.[ 6 ]for 
case D. 

Self-absorption within the detached shell has a limited 
effect, with a reduction ranging from 1% (model A) to 20% 
(model D), as compared to the optically thin approxima¬ 
tion (Fig. |4j. However, taking into account the background 
introduces a much larger effect (Fig. 0. 

The results depend on the adopted parameters in 
the model (mainly internal radius, expansion velocity and 
age). Smaller internal radius and expansion velocity, and/or 
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Figure 4. Hi line profiles of single detached shells for vari¬ 
ous circumstellar masses (A: 0.05 Mq, B: 0.1 Mq, C: 0.2 Mg, D: 
0.4 Mq,), no background. 



Velocity (km/s) 


Figure 5. Hi line profiles of a single detached shell 
(scenario 2, case D), and for various background levels 
{T bg = 0 K, 10 K, 30 K, 50 K). 

longer age would lower the average temperature in the 
detached shell. This would increase the effect of self¬ 
absorption, as well as that of the background absorption. 
The line profiles simulated with the A and B-cases repre¬ 
sented in Fig. U provide a good approximation to several 
observed HI line profiles (Libert et al. 2007, 2010, Matthews 
et al. 2013). 

As an illustration, we reproduce on Fig. [7] the spatially 
integrated profile of Y CVn observed by Libert et al. (2007) 
together with a recent fit obtained by Hoai (2015). For 


Figure 6. Density, velocity, and temperature profiles for a de¬ 
tached shell model (scenario 2, case D). 


Table 1. Model parameters (scenario 2), d = 200 pc, V eX p = 
8kms _1 , and M = 10 -7 Mq yr -1 . 


case 

age (yr) 

r f (arcmin) 

T/ (K) 

MflT.CS (M©) 

A 

5 x10 s 

3.85 

135 

0.05 

B 

10® 

4.14 

87 

0.1 

C 

2xl0 6 

4.47 

55 

0.2 

D 

4x10® 

4.83 

35 

0.4 


this fit, a distance of 321 pc (van Leeuwen 2007), a mass 
loss rate of 1.3xl0 -7 Mgyr -1 , and a duration of 7xl0 5 
years have been adopted. These parameters differ from those 
adopted by Matthews et al. (2013), who assumed 1.7xl0~ 7 
Mq yr -1 and a distance of 272 pc (Knapp et al. 2003). How¬ 
ever, by adopting a lower mass loss rate, and conversely, a 
longer duration, Hoai (2015) can fit the spatially resolved 
spectra obtained by the VLA and solve the problem faced 
by Matthews et al. at small radii. A difference between the 
mass loss rate estimated from CO observations and that 
adopted in the model may have several reasons, for instance 
an inadequate CO/H abundance ratio. 

4.3 Villaver et al. model (scenario 3) 

Villaver et al. (2002) have modeled the dynamical evolu¬ 
tion of circumstellar shells around AGB stars. The tempo¬ 
ral variations of the stellar winds are taken from the stellar 
evolutionary models of Vassiliadis & Wood (1993). 

For our Hi simulations, we used the 1.5-Mq circum¬ 
stellar shell models of Villaver et al. (2002) at various times 
of the thermally pulsing AGB (TP-AGB) evolution. We se¬ 
lected the epochs at 5.0, 6.5, and 8.0xl0 5 years, which corre¬ 
spond to the first two thermal pulses, and then to the end of 
the fifth (and last) thermal pulse. The density, velocity and 
temperature profiles are illustrated in Fig.[8] For these mod- 
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Figure 7. Y CVn integrated spectrum (Libert et al. 2007) 
and fit obtained by Hoai (2015) with scenario 2 (d = 321pc, 
M = 1.3xl0 -7 Mgyr -1 , age = 7xl0 5 years). 



Figure 9. HI line profiles of a circumstellar shell model around a 
1.5-Mq star during the evolution on the TP-AGB phase (5.0, 6.5, 
8.0x10 s years, Villaver et al. 2002), no background. The distance 
is set at 1000 pc. The first two profiles have been scaled by 37.7 
and 3.87, respectively, in order to help the comparison between 
the different line profiles. 


els, which can reach a large size (with radii of 0.75, 1.66 and 
2.5 pc, respectively), we adopt a distance of 1 kpc (implying 
a diameter of up to 17 / ). 

The results are shown in Fig. [9] In this scenario, the 
temperature in the circumstellar environment is maintained 
at high values due to the interactions between the successive 
shells (Fig. [Sj. The shape of the line profile is thus domi¬ 


nated by thermal broadening, and does not depend much 
on the epoch which is considered (although the intensity of 
the emission depends strongly on time, together with the 
quantity of matter expelled by the star). 

For the same reason, these results do not depend much 
on the background (<5% for Tbg=100K). Indeed in the 
models the temperature of the gas in the circumstellar shell 
always stays at a high level (> 10 3 K, except close to the 
central star in the freely expanding region). 

The predictions obtained with this scenario, in which 
wind-wind interactions are taken into account, differ clearly 
from those obtained with the previous scenario, in which 
the detached shell is assumed to result from a long-duration 
stationary process, by a much larger width of the line 
profiles (full-width at half-maximum, FWHM~ 16kms -1 ). 
This large width in the simulations for scenario n°3 results 
mainly from the thermal broadening, and also, but to a lesser 
extent, from the kinematic broadening (cf. Fig. [5J. 


5 DISCUSSION 

5.1 Optically thin approximation 

If absorption can be neglected, the intensity becomes pro¬ 
portional to the column density of hydrogen. For a source 
at a distance d, the mass in atomic hydrogen (Mm) can be 
derived from the integrated flux density through the stan¬ 
dard relation (e.g. Knapp & Bowers 1983): 

Mm = 2.36 x 10 7 d 2 J SmdV, 

in which d is expressed in pc, V in kms -1 , Shi in Jy, and 
Mm in solar masses (Mg). 

Our calculations allow us to estimate the error in the 
derived Hi mass of circumstellar envelopes that is incurred 
from the assumption that the emission is optically thin and 
not affected by the background. As an example, in Fig. HOI 
we show the ratio between the estimated mass (using the 
standard relation) and the exact mass in atomic hydrogen. 
The case without background illustrates the effect of self¬ 
absorption within the circumstellar shell for different mass 
loss rates. We adopt a freely expanding wind with a temper¬ 
ature profile in r -0,7 (as in Sect. 14.111 . with r expressed in 
arcmin., or a constant temperature (5K, 10 K, 20 K). The 
ratio clearly decreases with decreasing temperature in the 
circumstellar shell, increasing mass loss rate and increasing 
background temperature. 

In the constant temperature case with T = 5 K and 
Tbg= 5 K, the line profiles should be flat (cf. radiative trans¬ 
fer equation in Sect. and thus the “estimated” masses, 
exactly null. Our numerical calculations agree with this pre¬ 
diction to better than 3xl0 -3 , for mass loss rates up to 
10 -4 Mq yr -1 . 

The standard relation used for estimating the mass in 
atomic hydrogen should obviously be handled with caution 
in the case of the freely expanding wind scenario (n° 1). On 
the other hand, our calculations show that the deviation 
is much smaller for the two other scenarios (and basically 
negligible for scenario n° 3). This is mainly an effect of the 
high temperature in the detached shells resulting from the 
wind-wind interactions. 
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Figure 8. Density, velocity, and temperature profiles for the Villaver et al. (2002) model (scenario 3) at three different epochs (5.0 xlO 5 
years (left), 6.5 xlO 5 years (center), 8.0 xlO 5 years (right). 
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Figure 10. Ratio between the “estimated” mass in atomic hy¬ 
drogen and the real mass for the freely expanding wind case 
(scenario n° 1) with mass loss rates ranging from 10 —7 to 10 —4 
M 0 yr —1 (see Sect. 15. II) and different cases of temperature depen¬ 
dence (see text). Upper panel: no background. Lower panel: with 
a 5-K background. 

5.2 Spectral variations of the background 

The Hi absorption produced by cold galactic gas in the 
foreground of bright background emission may be shifted 
towards the velocity with highest background (cf. Levin¬ 
son & Brown 1980). To illustrate this effect in the case of 
circumstellar shells, in Fig. 1111 we show the results of our 
simulations for a 10~ 5 MQyr _1 freely expanding wind, as 
in Sect. 14. ll and a background temperature varying linearly 
between 10 K at —10 kms -1 , and 5K at +10 kms -1 . The 
absorption is clearly shifted towards velocities with high¬ 
est background. One notes also that the emission is shifted 
towards velocities with lowest background. 


Figure 11. Effect of a background intensity varying linearly 
from 10 K to 5 K across the line profile for a scenario 1 model 
with M=10 -5 Mq yr -1 . The curves labelled 'TgQ = 5K’, and 
’TgQ= 10K’, are reproduced from Fig. [3] 

In case of a intense and spectrally structured back¬ 
ground, some care should be exercised when comparing the 
Hi line centroids with the velocities determined from other 
lines. 


5.3 Comparison with observations 

Freely expanding winds have been definitively detected in 
the HI line in only two red giants: Y CVn (Le Bertre 
& Gerard 2004) and Betelgeuse (Bowers & Knapp 1987). 
The corresponding emission is relatively weak and diffi¬ 
cult to detect. Data obtained at high spatial resolution re¬ 
veal a double-horn profile (e.g. Bowers & Knapp 1987). 
It is worth noting that a high velocity expanding wind 
(Vexp ~ 35kms -1 ) has also been detected around the classi- 
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cal Cepheid S Cep (Matthews et al. 2012). A pedestal is sus¬ 
pected in a few HI line profiles that could be due the freely 
expanding region (Gerard & Le Bertre 2006, Matthews et 
al. 2013). The first scenario might also be interesting for in¬ 
terpreting sources in their early phase of mass loss, or for 
sources, at large distance from the Galactic Plane, embed¬ 
ded in a low pressure ISM. 

In general, sources which, up to now, have been de¬ 
tected in HI show quasi-gaussian line profiles of FWHM ~ 2- 
5kms _1 (Gerard & Le Bertre 2006, Matthews et al. 2013), 
a property which reveals the presence of slowed-down de¬ 
tached shells. These profiles are well reproduced by simula¬ 
tions based on the scenario n°2 presented in Sect. 14.21 assum¬ 
ing mass loss rates of a few 10 ' Mq yr - , and durations of 
a few 10 5 years. In particular, for Y CVn and Betelgeuse, the 
main Hi component has a narrow line profile (~3kms _1 ) 
and is well reproduced by this kind of simulation (Libert et 
al. 2007, Le Bertre et al. 2012). 

Sources with large mass loss rates (^ 5 x 1CU 7 Mq yr -1 ) 
have rarely been detected (with the notable exceptions of 
IRC+10216 and AFGL3068, see below). The simulations 
presented in Sect. [+3] show the line profiles that sources, 
such as those predicted by Villaver et al. (2002), should 
exhibit at the end of the thermal-pulse phase, with large 
mass loss rates, and with interaction with the local ISM. In 
these models, in which the evolution of the central star is 
integrated, the circumstellar envelopes result from several 
interacting shells, as well as from the ISM matter which has 
been swept-up. Shocks between successive shells maintain a 
high gas temperature (~4000K). 

For these models the calculated line profiles are not seri¬ 
ously affected by the background level, and the flux densities 
are large enough for allowing a detection up to a few kpc. For 
instance, in the GALFA-Hl survey (Peek et al. 2011), the 
3-cr detection limit for a point source in a lkms^ 1 channel 
is ~ 30 mjy. Saul et al. (2012) have detected many compact 
isolated sources in this survey. However, at the present stage, 
none could be associated with an evolved star (Begum et al. 
2010 ). 

Furthermore, several sources with high mass loss rates, 
such as IRC+10011 (WX Psc), IKTau (NMLTau) or 
AFGL 3099 (IZ Peg) which are observed at high galactic 
latitude, with an expected low interstellar Hi background, 
remain undetected (Gerard & Le Bertre 2006, Matthews et 
al. 2013). The simulations that we have performed based on 
the three different scenarios considered in this work cannot 
account for such a result. It seems therefore that, in sources 
with large mass loss rates (^ 5 x 10~ 7 MQyr _1 ), hydrogen 
is generally not in atomic, but rather in molecular form. 

Glassgold & Huggins (1983) have discussed the H/H 2 
ratio in the atmospheres of red giants. They find that for 
stars with photospheric temperature T*^ 2500 K, most of 
the hydrogen should be in atomic form, and the reverse for 
T*^ 2500 K. Winters et al. (2000) find that there is an 
anti-correlation between T* and the mass loss experienced 
by long period variables. It seems likely that stars having 
a mass loss rate larger than a few lO^Mgyr^ 1 have also 
generally a low photospheric temperature, with T* ^ 2500 K, 
and thus a wind in which hydrogen is mostly molecular. 

Recently, Matthews et al. (2015) have reported the de¬ 
tection of atomic hydrogen in the circumstellar environment 
of IRC+10216, a prototype of a mass losing AGB star at 


the end of its evolution with M~ 2xl0 -5 Mgyr -1 . The 
observed morphology, with a complete ring of emission, is 
in agreement with the predictions of Villaver et al. (2002, 
2012). They find that atomic hydrogen represents only a 
small fraction of the expected mass of the circumstellar en¬ 
vironment (< 1 %), supporting a composition dominated by 
molecular hydrogen. Unfortunately, a reliable line profile 
could not be extracted due to the low level of the emission 
and to a patchy background. The detection of Hi over a 
spectral range ~ lOkms^ 1 suggests a line width larger than 
commonly observed in evolved stars, which would make it 
compatible with scenario n°3. Gerard & Le Bertre (2006) 
have reported the possible detection of AFGL 3068, another 
carbon star with high mass loss rate (~ 10 -4 Mq yr -1 ). In 
this case, also, the line width (~30kms _1 ) is larger than 
expected for scenario n°2, and might be better explained by 
scenario 11 ° 3. Another possibility for this source which is at 
a large distance from the Galactic Plane (z ~ 740 pc) would 
be that we are mostly detecting a freely expanding wind not 
slowed down by its local ISM (i.e. scenario n°l). 

If the atomic hydrogen is of atmospheric origin (a frac¬ 
tion of 1 % is expected for a star with an effective temper¬ 
ature of 2200 K, Glassgold & Huggins 1983), its abundance 
should correspondingly be scaled down in our radiative 
transfer simulations. The effect of the optical depth on the 
line profiles could be considerably reduced for such a case. 
For stars with lower effective temperature (T*^2200K), 
atomic hydrogen might also be present in the external re¬ 
gions of circumstellar envelopes as a result of the photodis¬ 
sociation of molecular hydrogen by UV photons from the 
interstellar radiation field (Morris & Jura 1983). 

5.4 Case of a resolved source 

We have concentrated our study on the prediction of spa¬ 
tially integrated spectra. However, circumstellar envelopes 
may reach a large size (~ 2-3 pc, Villaver et al. 2002), and 
thus have a large extent over the sky. Also, interferometers 
provide a larger spatial resolution than single-dish antennas. 
It is thus interesting to examine how the line profile may vary 
as a function of position. As an example, in Fig. 1121 we show 
a spectral map that would be obtained for a detached shell 
observed over a background with Tbg— 50 K (scenario 2). 
The line appears mostly in emission and, as expected, de¬ 
lineates the detached shell. However, in the case of a high 
background level, the line appears also in absorption, in par¬ 
ticular in the external part of the detached shell where the 
lines of sight cross regions with gas at low temperature. 

Spatially resolved H 1 studies, with a careful subtrac¬ 
tion of the background emission, may thus reveal spectral 
signatures that hold information on their physical condi¬ 
tions. Such signatures could help to constrain the physical 
properties of the gas in a region where molecules are absent 
or not detectable. 


6 PROSPECTS 

We have simulated Hi 21-cm line profiles for mass losing 
AGB stars expected for different scenarios assuming spher¬ 
ical symmetry. 

However, AGB sources are moving through the ISM 
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Figure 12. Hi spectral map for a detached shell (case D in Ta- 
blc[l]with Tbg= 50 K), assuming a gaussian beam of FWHM = 1'. 
Steps are 1' in both directions. 

and their shells may be partially stripped by ram pressure 
(Villaver et al. 2003, 2012). As a consequence of the interac¬ 
tion a bow-shock shape appears in the direction of the move¬ 
ment, but also a cometary tail is formed which is fed directly 
from the stellar wind and from material stripped away from 
the bow shock. The cooling function and the temperature 
assumed for the wind have an important effect on the for¬ 
mation of the tail as shown in Villaver et al. (2012). Higher 
density regions formed behind the star will cool more effi¬ 
ciently and will collapse against the ISM pressure, allowing 
the formation of narrow tails. 

Gerard & Le Bertre (2006) have reported shifts of the 
Hi emission in velocity as well as in position for several 
sources. Matthews et al. (2008) have reported a shift in ve¬ 
locity for different positions along the tail of Mira (see also 
X Her, Matthews et al. 2011). These effects can also affect 
the HI line profiles, and thus the detectability. In addition, 
material lost by the AGB star should be spread along a tail 
that may reach a length of 4 pc, as in the exceptional case 
of Mira. On the other hand, Villaver et al. (2012) show that 
for sources with large mass loss rates at the end of their evo¬ 
lution, dense shells could still be found close to the present 
star position. 

We have assumed a background with a constant bright¬ 
ness. Of course, as explained in Sect. E3 this applies only 
to the cosmic background, and to a lesser extent to the 
galactic continuum emission. It does not apply to the galac¬ 
tic Hi emission which may show spatial structures of vari¬ 
ous kinds. The resulting effect may be more complex than 
that simulated in Sect. [4] For instance, an absorption line 
could form preferentially at the position of a peak of galac¬ 
tic Hi emission (a radiative transfer effect). Such a phe¬ 
nomenon may affect the predictions presented in Sect. 15.41 


Therefore a good description of the background will also be 
needed to model the observed line profiles. Such input may 
be obtained through frequency-switched observations for the 
galactic Hi component, and through the surveys of the con¬ 
tinuum at 21 cm which are already available (see Sect. [3]). 

It has to be noted that, in the position-switched mode 
of observation, the intrinsic line profile of the stellar source 
can also be spoiled by the patchiness of the galactic back¬ 
ground emission (observational artifact). The main source of 
confusion is the galactic HI emission which is structured spa¬ 
tially and spectrally. The classical position-switched mode 
of observation is not always efficient to correct the 21-cm 
spectra from the emission that is not directly associated 
with the star. More sophisticated methods with 2D mapping 
might be needed for subtracting the contaminating emission 
in these cases. Interferometric observations have the advan¬ 
tage of filtering the large-scale galactic emission. However, 
one should care that an intrinsic circumstellar emission is 
not also subtracted in this mode of observation. Also some 
artifacts may arise from incomplete spatial sampling of the 
large scale emission, as illustrated by the case of TX Psc 
(Matthews et al. 2013). If feasible, an excellent u-v coverage 
combined with maps from a single dish telescope providing 
small spacings has to be obtained. Also, for circumstellar 
shells angularly larger than the primary beam of the inter¬ 
ferometer, mosaicked observations combined with single dish 
maps are needed. 

Another caveat is that, when the distance to the source 
becomes larger, the foreground ISM material may play the 
role of an absorbing layer of growing importance (Zuckerman 
et al. 1980). The circumstellar shell line profile may thus be 
distorted by absorption due to the foreground cold material 
that shares the same radial velocity range. 

Sources with high mass loss rate (~ 10 -6 - 
10 _5 Mgyr _1 ) tend to be concentrated towards the 
galactic plane. They are expected to dominate the con¬ 
tribution of AGB stars to the replenishment of the ISM 
(LeBertre et al. 2003). The recent detection of IRC+10216 
by Matthews et al. (2015) shows that atomic hydrogen 
should be present in these sources and that the Hi line 
at 21 cm can be used to probe the morphology and the 
kinematics of stellar matter decelerated at large distance 
from the central star. However, as discussed above, when 
the background is large, a proper modelling of the line 
profiles will be necessary. 


7 SUMMARY AND CONCLUSIONS 

We have simulated HI 21-cm line profiles expected for sev¬ 
eral different scenarios representing different evolutionary 
stages of evolved stars, and thus corresponding to different 
AGB circumstellar structures. We have relaxed the optically 
thin hypothesis which was assumed in previous works, and 
included the emission from the background. 

Self-absorption may be important in freely expanding 
circumstellar shells, as well as in some detached shells re¬ 
sulting from the interaction of the stellar winds with the 
local ISM. The Hi line profile may also be affected by the 
background level and by the spectral profile of this back¬ 
ground emission. 

The numerical simulations that we have performed show 
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that, under certain conditions, the observed Hi 21-cm flux 
densities from mass-losing stars can be significantly reduced 
by taking into account optical depth effects and the presence 
of the background emission, but not to such a level such as 
to account for the non-detection of several sources. There¬ 
fore, one should consider that molecular hydrogen instead 
of atomic hydrogen likely dominates in sources with high 
mass loss rates few 10 -7 Mq yr _1 ), probably an effect of 
their low atmospheric temperature. Still, the recent results 
of Matthews et al. (2015) show that the Hi line at 21 cm 
can be a useful probe of the outer regions of sources with 
low stellar effective temperature (<2500K). 

For sources with mass loss rates ~ 10 _ ' Mq yr _1 , which 
are detected in Hi, the global agreement between the ob¬ 
served line profiles and the simulations based on the second 
scenario suggests that their central stars undergo mass loss 
smoothly over several 10 5 years. 
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